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We present first measurements of the (j>- meson elliptic flow (v2(pt)) and high statistics pr distri- 
butions for different centralities from ^/sjvjv = 200 GeV Au+Au collisions at RHIC. In minimum 
bias collisions the V2 of the cj> meson is consistent with the trend observed for mesons. The ratio 
of the yields of the Q to those of the as a function of transverse momentum is consistent with 
a model based on the recombination of thermal s quarks up to pt ~ 4 GeV/c, but disagrees at 
higher momenta. The nuclear modification factor (Rcp) of cf> follows the trend observed in the K% 
mesons rather than in A baryons, supporting baryon-meson scaling. Since </>-mesons are made via 
coalescence of seemingly thermalized s quarks in central Au+Au collisions, the observations imply 
hot and dense matter with partonic collectivity has been formed at RHIC. 
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The primary aim of ultra-rclativistic heavy-ion colli- 
sions is to produce and study a state of high-density 
nuclear matter called the Quark-Gluon Plasma (QGP), 
the existence of which is supported by lattice QCD cal- 
culations [J 0, 9- In the search for this new form 
of matter, penetrating probes are essential in order to 
gain information from the earliest stage of the collisions. 
Phenomenological analysis [3| has suggested a relatively 
small hadronic interaction cross section for ^-mesons al- 
though discussions about the <p— proton interaction cross 
section are yet to be conclusive [JIB] ■ Therefore 0-mesons 
from high-energy nuclear collisions are expected to pro- 
vide information about the early partonic stages of the 
system's evolution since they should remain mostly unaf- 
fected by hadronic interactions. This is further supported 
by recent measurements Q which have ruled out the idea 
of </>-meson production by kaon coalescence. 

Elliptic flow, w 2 , is an observable which is thought to 
reflect conditions from the early stage of the collision 
[1, Q. In non-central heavy-ion collisions, the initial spa- 
tial anisotropy of the overlap region of the colliding nuclei 
is transformed into an anisotropy in momentum space 
through interactions between the particles. Systematic 
measurements of the V2 for the strange hadrons Kg, A, 
S and f2 suggest that collectivity is developed at the par- 
tonic stage at RHIC [T3, EJ ■ The evidence for partonic 
collectivity, one of the conditions for QGP formation, will 
be further strengthened if it can be shown that 0-mesons 
flow like the other mesons. 

A mass ordering predicted by hydrodynamics [HI, [HI, 
[Til l for i>2 (pt) of identified particles has been observed 
for px < 2 GeV/c. At intermediate transverse momen- 
tum, 2 < pt < 5 GeV/c, a separation of baryons and 
mesons has been observed in measurements of both v-i 
and the nuclear modification factor, Rcp [ill . [H, 
These results are consistent with calculations from quark 
recombination models [Tt], EH, El, H3] implying the de- 
confinemcnt of the system prior to hadronization. The 
cj> is a vector meson, comparable in mass to the proton 
and A baryons with a relatively long lifetime. Its V2 and 
Rcp will provide a critical test of the assumed underly- 
ing dynamics. In addition, as argued in [2l| . the ratio 
of the fi-baryon over 0-mcson yields can be used to test 
the nature of light-quark thermalization in the medium. 
The model predicts that the ratio of the Q to <j) yields 
will rise monotonically. 

The results presented in this paper were obtained with 
the STAR detector [22| at the Rclativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory. 
The detector components used in this analysis were the 
Time Projection Chamber (TPC), and trigger detectors, 
namely the zero degree calorimeters. Central collisions 
were selected using the Central Trigger Barrel. The col- 
lision centrality was determined by the charged hadron 
multiplicity within pseudo-rapidity |t^| < 0.5. 

High-statistics Au+Au data were taken in the 2004 
run at ^Jsnn = 200 GeV. The minimum bias dataset 
used in this analysis consisted of ^13.5 million events and 



the central-triggered dataset comprised about 10 million 
events. Events were required to have a primary vertex 
z position (where z is the direction of the beam axis) 
within 30 cm of the center of the TPC. Events from the 
minimum bias dataset were divided into 8 centrality bins: 
0-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60- 
70%, and 70-80% of the measured cross-section. The 
central-triggered dataset was used to extract the 0-5% 
and 0-12% data. 

The 4> yield in each pt bin was extracted from the 
invariant mass (rrii nv ) distributions of K + + K~ candi- 
dates after subtraction of combinatorial background esti- 
mated using event mixing Q . The kaons were identified 
through their dE/dx energy loss in the STAR TPC [II]. 
Including the detector resolution, the values of the recon- 
structed^ mass and width are consistent with the PDG 
values [23j. The relative systematic uncertainty due to 
the dE/dx cut was estimated to be ~ 8% by using differ- 
ent cuts and comparing the yields after a particle identi- 
fication efficiency correction. Uncertainty in the residual 
background shape of the m m „ distributions resulted in 
a contribution of about 4.5% to the errors on the final 
yields. 

The 0-meson vi results were obtained using the V2 vs. 
rriinv method described in ref. [24j . The method involves 
calculating the vi of the same-event distribution as a 
function of rrii nv and then fitting the resulting V2(rrii nv ) 
distribution using: 

v 2 {rriinv ) = V2sa{m inv ) + v 2 b (iTlinv 

)[i )] (i) 

where V2s = V24, is the signal V2 and V2B is the back- 
ground V2- a{rrii nv ) = S/(S + B) is the ratio of the signal 
over the sum of the signal plus background of the m, nnj 
distributions. It was extracted from fits (Breit-Wigner 
plus a linear function) to the <f> mass-peak for each pt 
bin. For each pt bin, the V2(rrii nv ) was fitted using Eq. [T] 
in order to extract the fitting parameter V2S and V2B 
was parameterized using a linear or quadratic function 
in rriinv These results are consistent with results using 
an established method [25[ where the 0-meson yield is 
plotted as a function of the difference between its az- 
imuthal angle and the estimated reaction plane angle, 
(<fi — ^f). The values of i>2 are extracted from the fitting 
to the function dN/dcj) = P (l + 2v 2 cos(2(0 - *))). 

In the top panel of Fig. Q] we present the first mea- 
surement of the differential elliptic flow, V2(j>t), of the 
0-meson from Au+Au collisions for four centrality bins. 
In this and the following figures, the vertical error bars 
on the 4> data points indicate the statistical errors while 
the shaded bands indicate the extent of the systematic 
uncertainties. The systematic errors vary from point to 
point including uncertainties in extracting the signal for 
obtaining a(rrii nv ) and differences in the reaction plane 
resolution determination. For minimum bias collisions, 
an additional contribution to account for the different 
methods of extracting the V2{pt) values is also included 
in the systematic error. Non-flow effects [26|, [27| arc not 
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FIG. 1: (color online) Top panel: The elliptic flow, V2(pr), 
for the 0-meson as a function of centrality. The vertical 
error bars represent the statistical errors while the shaded 
bands represent the systematic uncertainties. For clarity, 
data points are shifted slightly. Bottom panel: Minimum 
bias V2(pt) for the <f>- meson compared to results for A and 
Ks EH] ■ The dashed and dotted lines represent parameteriza- 
tions inspired by number-of-quark scaling ideas from ref. 
for NQ=2 and NQ=3 respectively. 



included in the systematic error. As expected, V2(pt) in- 
creases with increasing eccentricity (decreasing central- 
ity) of the initial overlap region. This trend is also illus- 
trated in Table U which presents the pr-integrated values 
of 0- meson elliptic flow, (V2), calculated by convoluting 
the V2(pt) with the respective px spectrum for three cen- 
trality bins. It should be noted that the centrality depen- 
dence of the (^2) of ^-mesons is consistent with that of 
charged hadrons [27j | . 



TABLE I: Integrated elliptic flow, (172), for the 0-meson for 
three centrality bins. 



Centrality (%) 



<ga) (%) 



40 - 80 
10 - 40 
0-5 



..5 ± 



l.l(stat) 
0.2<sys) 

0.2(sys) 
2 -j^ 1 1.2(stat) 
O.b(sys) 



The lower panel of Fig. [TJ shows the minimum 
bias (0-80%) result compared to parametcrizations for 
number-of-quark scaling for mesons (NQ=2) and baryons 



(NQ=3) whose free parameters have been fixed by fit- 
ting to the A and K s results simultaneously [28[ . In this 
case, for pt < 2 GeV/c, the cf> V2 follows a mass-ordered 
hierarchy where the values of V2, within errors, fall be- 
tween those of the heavier A (open circles) and lighter 
Kg (open-squares). However, at intermediate pt, be- 
tween 2-5 GeV/c, the <j> V2 appears to follow the same 
trend as K s . When we fit the i>2_(pt) of mesons with 
the quark number scaling ansatz [28j , the resulting fit pa- 
rameter NQ = 2.3±0.4. The fact that the </> V2(pt) is the 
same as that of other mesons indicates that the heavier s 
quarks flow as strongly as the lighter u and d quarks. As 
previously mentioned, (^-mesons are not formed through 
kaon coalescence Q and do not participate strongly in 
hadronic interactions. Therefore the results demonstrate 
partonic collectivity 

Figure [2] shows the px distributions of ^-mesons as 
a function of centrality The central-triggered dataset 
was used to obtain the most central spectrum while the 
other distributions were obtained using the minimum 
bias dataset. The error bars shown in Fig. [2] are sta- 
tistical only. In the figure, the errors are smaller than 
the size of the data points. 
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FIG. 2: (color online) Transverse momentum distributions of 
</>-mesons from Au+Au collisions at ^/snn = 200 GeV. For 
clarity, distributions for different centralities are scaled by 
factors of ten. Dashed lines represent the exponential fits to 
the distributions and the dotted lines are Levy function fits. 
Error bars represent statistical errors only. 

Each pt spectrum in Fig. [2] has been fitted using 
both an exponential function (dashed lines) in and 
a Levy function (dotted lines) which has an exponential- 
like shape at low px and is power-law-like at higher px- 
While the central data are fitted equally well by both 
functions the more peripheral spectra are better fitted 
by the Levy function indicating less thermal contribu- 
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tions in peripheral collisions. 
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FIG. 3: (color online) The N(Q,)/N(<f>) ratio vs. p T for three 
centrality bins in ^/sjvjv = 200 GeV Au+Au collisions. The 
solid and dashed lines represent recombination model predic- 
tions for central collisions [2l[ for total and thermal contribu- 
tions, respectively. 

In Fig. [3l the ratios of N(fl) /N(<p) vs. pp are presented 
as a function of centrality. The fl datapoints are from 
ref. HI (for 0-10%) and ref. [13] for the other central- 
izes. The errors of the ratios are dominated by the 
datapoints. Also shown in the figure are recombination 
model expectations for central collisions [2l[ based on <fr 
and production from coalescence of thermal s quarks in 
the medium. The model describes the trend of the data 
up to pt ~ 4 GeV/c but fails at higher p?. Other mod- 
els based on dynamical recombination of quark s fl9l . l3l| 
were also compared to the data. However, ref. over- 
predicts the ratio while ref. [3l| gives the wrong shape. 
With decreasing centrality, the observed N(n)/N((f>) ra- 
tios seem to turn over at successively lower values of pt 
indicating a smaller contribution from thermal quark co- 
alescence in more peripheral collisions. This is also re- 
flected in the smooth evolution of the spectra shapes from 
the thermal-like exponential to power-law shapes shown 
in Fig. |H 

The nuclear modification factor Rcp(pt) measures the 
change of px distributions from peripheral to central col- 
lisions and has been measured for most of the identified 
hadrons. In Fig.dJ the high statistics </>-meson Rcp (solid 
circles) is compared to Kg (open triangles) and A (open 
squares) from ref. [ll[ for two different centrality com- 
binations (upper and lower panels). In both panels the 
binary-scaled yield of 0- mesons is suppressed (Rc p below 
unity) in central compared to peripheral collisions. It has 
been shown that u, d and s quarks are approaching equili- 
bration at hadronization [291 . The <b -meson Rc p is more 
consistent with that of K s (meson) than of A (baryon) for 
the 0-5%/40-60% case (upper panel). This is as predicted 
by particle production models based on recombination of 
thermal quarks fljil ]. For the more peripheral bin (lower 
panel) the <fi Rcp falls between that of the A and K s . 
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FIG. 4: (color online) The Rcp of mid-rapidity <^>- mesons 
produced in ^/snn = 200 GeV Au+Au collisions: (top) 0-5% 
vs. 40-60% and (bottom) 0-5% vs. 60-80% The shaded bands 
represent the uncertainties in the Glauber model calculations 
for (N bin ) and {N par t) [12] ■ Also shown are results for A and 
Ks El and protons and tt + 



In the 60-80% centrality bin (see lower, binary collision- 
scaled (f> production is very similar to that in p+p and 
d+Au collisions where strangeness production is canoni- 
cally suppressed (34[. Therefore a baryon- meson scaling 
behaviour of Rcp is not expected in the lower panel of 
Fig. 2J In addition, for baryons and mesons respectively, 
there seems to be an ordering in terms of strangeness 
content. This has also been observed in Raa for strange 
particles [35| . 

In summary, we have presented first measurements of 
the elliptic flow of </>-mesons as a function of collision cen- 
trality in Au+Au collisions at ^/snn = 200 GeV. At low 
Pt (< 2 GeV/c), V2 is consistent with hydrodynamical 
expectations. At intermediate px (2 < pp < 5 GeV/c), 
V2 of 0-mcsons is consistent with numbcr-of-quark scal- 
ing for mesons. These observations indicate the develop- 
ment of partonic collectivity in the medium. Measure- 
ments of the <j> Pt spectra as a function of centrality 
show an evolution of the spectral shape from exponen- 
tial to power-law-like with decreasing centrality, reflect- 
ing the increasing contributions from hard and possibly 
other non-equilibrium processes in more peripheral colli- 
sions. The result of a recombination model [2l| is consis- 
tent with the trend of the central N({i)/N(<f>) ratio up to 
Pt ~ 4 GeV/c which covers more than 95% of the hadron 
yields. At higher px, the model fails. The </>-meson Rcp 
resembles the K s for the 0-5%/40-60% case which is con- 
sistent with meson scaling. Since 0-mesons are made 
via coalescence of seemingly thermalized s quarks in cen- 
tral Au+Au collisions, the observations imply hot and 
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dense matter with partonic collectivity has been formed 
at RHIC. 
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